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Abstract 

In 1991, Ursell gave a strong form of Watson’s lemma for the Laplace integral 

e~^*f{t) dt (x —>■ +oo) 

in which the amplitude function f(t) is regular at the origin and possesses a Maclaurin 
expansion valid in 0 < t < i?. He showed that if the asymptotic series for the 
integral as a: —> +oo is truncated after rx terms, where 0 < r < R, then the resulting 
remainder is exponentially small of order In this note we extend this result 

to include situations when f{t) has a branch point at t = 0 and when x is a complex 
variable satisfying | arg x| < 7r/2. 
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1. Introduction 

Suppose that the function F{z) possesses the asymptotic expansion 

OO 

( 1 . 1 ) 

k=0 

as z ^ OO in a sector S of the complex z-plane. If we denote by Rn{z) the error in 
approximating F{z) by the sum of the first n terms, then we have exactly 

n—1 

F{z) = Y,CkZ-^+ Rn{z) (zES). 

/c=0 

By Poincare’s definition of an asymptotic series, the error (or remainder) term satisfies 
the bound \Rn{z)\ < An\z\~'^, where An is an assignable constant for z in the sector 
S. The integer n is arbitrary, but fixed: truncation of the series (1.1) at the nth term 
consequently yields an error that decays algebraically as z ^ oo in S. When \z\ is large, 
the successive terms in a typical asymptotic expansion (1.1) initially start to decrease 
in absolute value, reach a minimum and thereafter increase without bound given the 
divergent character of the full expansion. If, however, the series is truncated just before 
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this minimum modulus term is reached, then this process is called optimal truncation 
and the finite series that results is the optimally truncated expansion. As we shall see in 
a specific case below, the resulting rate of decay of Rn{z) is then greatly enhanced. 

To illustrate we consider the exponential integral 

roo p—xt 

£{x) := Ei{x) = dt, 

where x > 0, which provides one of the simplest examples of an asymptotic expansion 
complete with error bound. Substitution of the identity 


1 

1 +1 


n—1 


k=0 


i-tr 

i+t 


(t /1), 


followed by termwise integration, yields 


^1 U] roo fTL 

f(a:) = Ux) = {-r - 


OO ^Up — xt 


+ t 


dt. 


( 1 . 2 ) 


It is readily seen that 

roo Y) I 

\Rn{x)\< = — , 

Jo x"- 

so that this expansion enjoys the property that the remainder is bounded in magnitude 
by the first neglected term and has the same sign. 

For a given value of x, the smallest term in absolute value of the series in (1.2) occurs 
when k = [xj (except when x is an integer, in which case there are two equally small 
terms corresponding to k = x — 1 and k = x). If we denote the optimal truncation index 
by N, then we have 

^-1 L, 

k^O ^ 

It is clear that RNix) is a discontinuous function of x, since N changes each time x 
passes through an integer value. Use of Stirling’s formula N\ (27r)V2e-^Ar^+V2 ^ 
approximate N\ for large ~ x then produces the estimate 


N\ 


\RNix)\ <Zn- (27r)2 


1 e 


-N]\fN+^ 


X 


N 


~ (27rx)2e 


as X —?• +00. This shows that at optimal truncation the remainder for £(x) is of order 
xi/ 2 g-a; ^ _I_QQ consequently that evaluation of £{x) by this scheme will result 
in an error that is exponentially small. This level of asymptotic approximation has been 
given the neologism superasymptotics in [1]; sometimes it is referred to as exponential 
improvement. 

As a rule, it is found that similar exponential improvement at optimal truncation can 
be achieved in asymptotic expansions of other functions, most notably the confluent hy¬ 
pergeometric functions which include many well-known special functions, such as Bessel 
and Airy functions. In these cases the remainder can be written explicitly in the form 
of an integral and its value estimated to establish its exponentially small nature. The 
idea of terminating an asymptotic expansion near its smallest term and estimating the 
remainder is well known in the asymptotic theory of converging factors; see, for example, 
[7, §7.34] and [3, Ch. 14]. 
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There are, however, few general results of this character. In 1958, Jeffreys [2] consid¬ 
ered the Laplace transform of (with ^ > —1, /(O) 7 ^ 0) which has the expansion 

for X —)• -|-oo 


f e ^H>^f{t)dt='^Ur + Rn{x), Ur = (1.3) 
Jo r.x 

where Rn{x) is the remainder after n terms. If the singularity of f{t) closest to the origin 
is situated at Re'‘°‘ (a / 0), then the least term in modulus in the above series occurs 
when n ~ Rx. When n is so chosen, he showed that the optimal remainder is given by 

^ - 1 {x^+ 00 ). 

In the most common case we have a = tt and the above result shows that the remainder 
is then approximately half the least term in the expansion. 

A significant advance was made by Ursell [ 6 ] who established the following result. 

Theorem 1. Let f{t) be analytic in the disc \t\ < R with the Maclaurin expansion 

00 

f{t) = ^ Cnt^ (0 < t < R). 

n=0 

Let r be a fixed positive quantity such that 0 < r < R and suppose that \f{t)\ < Ke^^ 
when r < t < 00 , where K and fi are positive constants. Then 

Iix)= e-^^f{t)dt = J2^^ + ^nfix), (1.4) 

“'O n=0 ^ 

where, if n* is chosen to equal [rxj, the remainder Rnfix) = 0{e~'^^) as x ^ -boo. 

Ursell’s proof relied on simple bounds for the incomplete gamma functions and was 
particularly elegant and straightforward; see also [5, p. 76] for a detailed account. We 
observe that the series appearing in (1.3) and (1.4), when the upper limit is extended 
to inhnity, is the standard asymptotic expansion of the integrals on the left-hand sides 
as X —)• -boo obtained by application of Watson’s lemma; see, for example, [4, p. 44]. 
Ursell’s result shows that if we take 1 -b [rxJ terms in the asymptotic series, then the 
resulting approximation to I{x) is exponentially accurate, with the error being 0 (e“™) 
as X ^ -boo. 

Our aim in this note is to extend the result in Theorem 1 to situations where the 
function f[t) possesses a branch point at the origin and to include sectors of the complex 
x-plane. We carry this out using a modification of the procedure adopted by Ursell. 


2. The modified expansion 

We consider the function I{z) defined by the Laplace integral 

fOO 

I{z)= e-^^f{t)dt, (2.1) 

Jo 

where 2 : = xe*^ is a large complex variable with x > 0 and phase 9 = arg 2 : satisfying 
1^1 < ^vr — J, J > 0. The amplitude function f{t) satisfies the following conditions: (i) 
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f{t) is holomorphic in the sector S: —ai + 5 < arg t < a 2 — 6; (ii) f{t) possesses the 
absolutely convergent expansion 

°° n+/3 

= ^ {\t\ < R), 

n=0 

where ^ > 0 and 3f?(/3) > 0; and (hi) \f{t)\ < {t G S), where A and a are positive 

quantities. Then we have 

Theorem 2. Let f{t) satisfy conditions (i) - (Hi) above where we suppose that ai, a 2 > 
^TT. Let r be a fixed positive number sueh that r < R and n* = [ixr\z\ + /r — 3f?(/3)J. Then 

'W= ( 2 . 2 ) 

n<n* \ r' / 


where the remainder Rn, (z) satisfies the bound 

as \z\ GO in the seetor \ arg z\ < — 6, 6 > 0. 

Proof. We rotate the integration path in (2.1) by —6 to hnd 



roo 

e-^V(t) dt = e-^^ / e-^^fire-^^) dr, 

Jo 


where the contribution from the arc at infinity vanishes by assumption (iii) when \z\ > 
acosecS [3, p. 107]. Then, for 0 < r < R, 


I{z) 


e"*® y Cn r e-^^Te-y^+d)lu-^dT + J 
n=0 do 


E- 

n=0 


Cn 

y{n+p)lfj. ^ 


n + /3 
T 


, rx + d 


(2.3) 


where 7 ( 0 , x) is the lower incomplete gamma function and 


fOO 

J = e-^^ / e-^^fire-^^) dr. 

Jr 

The sum in (2.3) converges since 


— XT_ 

e T 


n+at(/3) 


-1 


dr < 


‘ njMSl-l 
T dr 


^j,(n+5R(/3))//i 
n + 3ft(/3) 


It may be remarked that (2.3) represents the first stage of a Hadamard expansion process 
as discussed in detail in the monograph [5, Ch. 2]. 

Let us now split the series in (2.3) at n = n*, where n* is to be specihed. Then upon 
use of the relation 

7(0, x) + r(a, x) = r(a) 

between the upper and lower incomplete gamma functions, we have 


I{z) 


SP i V f ^ ^ 

hi z^^^d)iu\ V 


-r 


n -\- j3 \ 1 Cn 

^ d d n>n. ^ 


n +13 
T 



+ J- 


(2.4) 
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If we now define the remainder (z) by 


Hz) = E 


n<n* 




n + P 


“I" -^n* 


it then follows from (2.4) that 


\RnM\ < E 


n<n* 

We now choose 


j(n+/3)/At| 


n + (3 


, rx 


+ E 

n>n* 


j(n+/3)/At| 


7 


n + (3 


, rx 


+ I'^l- 


= [lJ,rx + fi-^if3)\, (2.5) 

where we recall that /x > 0 and 3ft(/I) > 0 by hypothesis, and make use of the bounds 
given in (A.l) and (A.2) 


n + /3 
n + /3 


, rx 


-,rx 


< 2e-™(rx)("+*(^))/^ (0 < n < n*) 

< e-"^(rx)(”+^(^))/^ (n > n*). 


Thus we find 

\RnAz)\ < 2Ae-^^ E + E + 


n<n* 

oo 


n>n* 


< 2 Ae-""Elc-l^ 

n=0 


(n+5R(/3))//i 


+ l-^l) 


( 2 . 6 ) 


where A = exp[09(/3)/;u]. By condition (ii) it follows that X) is convergent 

for r < R. 

Let a > O' be a value of 3ft(z) for which /q°° e~^^f(t) dt converges. Then when x > a, 
we find by partial integration [3, p. 72] 

POO TOO 

J= e-^^f{Te-^%dT= / e-(’’’-“)"e-“V(re-*'^)dT 

Jr J r 

POO 

= [x — a) / F{t] 9 ) dr, 

Jr 

where 

F{T-e)= /E-“^/(Te-*®) dr. 

Jr 

If L{6) denotes the supremum of |T'(t; 0)| along the ray [r, oo) situated in S then^ 

POO 

I J| < (x - a)L{e) / = L(0) . (2.7) 

Jr 

By combining (2.6) and (2.7), we see that with the above choice of n* the remainder 
satisfies |ift„^(z)| = 0{e~'’^), thereby establishing the theorem. □ 

We now suppose that a singularity of f{t) lies on the circle of convergence jtj = R 
situated at to = Re~^^, where 0 < 'ip < An analogous treatment applies when 
there is a singularity at Re^^. During the process of path rotation, it is now possible to 
cross over the singularity, thereby receiving either a pole contribution or a branch point 
contribution depending on the nature of the singularity at to- 


'^A bound on | J| can also be obtained from condition (iii); see [6], [5, p. 76]. 
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Theorem 3. In the case when there is a singularity of f{t) on the circle of convergence 
|t| = R at the point to = Re~'^'^, where 0 < i/' < ^vr, then the remainder Rn,{z) satisfies 
the bound 

(z) = + T(0) (2.8) 

where n* is defined in (2.5) when r < R and 

fO —l7r + 5<9<fi 
T{e) = 2 “ “ r 

V' + (5<0<^7r — 5. 

An analogous result holds if the singularity is situated in the first quadrant at to = Re^'^. 


Proof. The treatment when 9 < if and when 9 > if follows that described in the 
proof of Theorem 2, where in the latter case there is the addition of the contribution 
from the singularity at to. In the case of a simple pole this yields a contribution of 
0 (exp {—zto)) = 0(exp (—|z|ii cos(0 — V’)))- In the case of a branch point singularity, the 
additional contribution is given by the loop integral 

/ e-^^f{to + u)du, 

J OO 


which is again of 0(exp (—zto))- 

When 9 = if,we can take the integration path for J to be the path commencing at 
and passing to infinity parallel to the real t-axis, thereby avoiding the singularity 
at to- Then 

J= e-^^f{t) dt = e-^^ / e-^“/(n + re”*’^) du, 

iV" JO 

where on the path u G [0, oo) we have by condition (hi) \ f{u + re“*^)| < ^ 

AeAu+r)_ Hence 


I J| < Ae 


roo 


Ae-rix-^l e- 


-rx I ^- 

0 »£(z) - a 


< R' 


X cos if ’ 


where A' is a constant. By the argument leading to (2.6), the remainder satisfies 
\RnSz)\ = 0(e“™) when 9 = if and n* is chosen according to (2.5), thereby estab¬ 
lishing the theorem. □ 

Remark 1. If there is a second singularity on |t| = R with | arg t| < ^vr then the modi¬ 
fication to Theorem 3 is straightforward. Similarly, if there is an additional singularity of 
f{f) beyond the circle of convergence \t\ = R situated at to = R'e~'’^, say, where R' > R 
and 0 < V’ < fTT, then the above argument is easily modified to yield the remainder 

{z)\ = 0(e-"l"l) + T(0) (2.9) 


Remark 2 . The second order estimate in (2.8) will only be significant (when 9 > if) if 

9 — if > aiccos{r/R). (2-10) 

A similar remark applies to (2.9) when R is replaced by R'. 
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3. Numerical verification 


We present some examples which illustrate the remainder estimates in Theorems 2 and 
3 when the underlying asymptotic series is truncated after n* terms, where n* is defined 
in (2.5). 


Example 1. We consider the confluent hypergeometric function 

U{a,a — b + l,z)= / (3ft(a) > 0) 

r(a) Jo 


for which 


fit) = ^ 1^1 ^ 


where {a)n = r(a + n)/r(a) is Pochhammer’s symbol. The function f{t) has a singularity 
at t = —1, so that ii = 1, and is associated with the parameters fj, = 1, /J = a. Then the 
remainder is dehned by 


Rn, (z) = U{a, a — b + 1, z) — z 


-“E 

n<nH 


{-na)n{b)r 


where the truncation index n* = rjzl with r < 1. According to Theorem 2, the remainder 
Rn*{z) = as \z\ —)■ oo in | arg z\<^ti — 5. 

It is known [4, §13.7(iii)] that when n* = \z\ the remainder for the U confluent 
hypergeometric function satisfies Rn,{z) = 0(e~^^^) uniformly for | arg z\ < vr. 


Example 2. The function 

where Ko( 2 ;) is the Struve function dehned in [4, Eq. (11.2.5)]. The function 
fit) = (1 + t^)-2 = ELEt t^n (1^1 < 

n 'n'- 
n=0 

has singularities at t = ±z and is associated with the parameters /U = /3 = ^. Then the 
remainder is dehned by 

R iz) - Hz) - - V 

where the truncation index n* = r| 2 ;| with r < 1. Erom Theorem 2, Rn*iz) = 0(e“^l^l) 
as l^l —)■ 00 in I arg z| < ^tt — 5. 

In Eig. 1 we show the variation of \Rn,iz)\ in Examples 1 and 2 as a function of 9 
in the range 0 < 0 < ^tt for different values of x. It can be seen that each curve lies 
close to the value of the order estimate and that (on the scale of the hgure) there 

is little variation in the value of \Rn,iz)\ over the range 0 < 0 < ^vr. Conjugate values 
are obtained for — ^vr < 6* < 0. 

Example 3. The function 

POO p fj't 

Hz) = l (o<^<H 
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Figure 1: Plots of \Rn, {z)\ against 0/tt (on a log^Q scale) for a; = 5,10,15, 20 and r = 0.8 for (a) 
the confluent hypergeometric function U{a,a — b + 1, z) with a = i, 6 = | and (b) the Struve 
function in Example 2. 


is associated with the amplitude function 

OO 

f{t) = e^^{l - (|t| < 1 ), 


n=0 


SO that ^ = (3 = 1. The function f{t) has a simple pole at to = e which produces the 
residue contribution 27ri exp (—Then, with the remainder defined by 

Rn,{z) = m- E n\{ze-^^)-^-\ 

n<n* 


we have from Theorem 3 that as \z\ —)• oo 


(3.1) 


Example 4. Our final example is the function 

r°° e~^^dt 
0 Vl - te*’/’ 


I{z) = [ 
do 


(0 < -0 < ^tt). 


Here we have 


OO /I 


f{t) = (1 - = E (|i| < 1 ), 

n\ 


n=0 


so that ^ = j3 = 1. The amplitude function f{t) possesses a square-root branch point at 
to = which produces the contribution 


dio-) e ^‘dt ^ 


Vl ~ te*d 7—00 (ue*d)2 V 

during the path rotation argument. With the remainder defined by 

RnAz) = I{z)- E 

n<n* 

we have the order estimate^ given in (3.1). 

^We omit the term in the second exponential. 
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Figure 2: Plots of logiQ(e''l^l |i?„, (z)]) against dj-K (solid curves) for r = 0.8 and two values of ij) 
when X = 20: (a) for Example 3 and (b) for Example 4. The dashed curves show the variation of 
log^o(6 '^'l^n. (^)l) s-s a function of 9/tt when i/i = O.Itt. 


In Fig. 2(a) we give an example of the variation of log;^Q(e^l^l |i?„^( 2 ;)|) for the integral 
in Example 3 as a function of 6 when x = 20 and tp = O.dvr. For this value of tp, 
the condition ( 2 . 10 ) is not satisfied and so the behaviour of \RnSz)\ is controlled by 
the 0(e“'’l^l) term. When pi = O.lvr, however, the condition (2.10) yields 9 — p > 
arccos 0.8 ~ 0.2057r. Thus, for 0 >O. 3 O 57 r the behaviour of \Rn,{z)\ is controlled by the 
pole contribution, as can be seen from the figure where the value changes by about five 
orders of magnitude. The dashed curve shows the variation of ghlcos( 6 »-' 0 ) it is 

seen that this curve approaches the value log^g 27r = 0.79818 confirming the dominance 
of the second term in (3.1) for 0 >O. 37 r. Fig. 2(b) shows the variation of the remainder 
in Example 4 for the same values of r, x and p as in Example 3. Similar remarks apply 
in this case. 


Appendix: Derivation of bonnds on the incomplete gamma functions 


Let io = a + ib, where a and b are real, and suppose that y > 0. Then, following Ursell 
[ 6 ], we have 

|r(u; + l,x)| = 


poo 

/ e-H‘^dt 

= g-X^«+l 

poo 

/ e-'^pi + uYdu 

Jx 


Jo 




e-^^il + uPdu. 


When — 1 < a < 0, we have 

poo 

|r(a; + 1, x)| < e->^x“+' / e-^^du < e^x“+' (x > !)■ 

Jo 

When 0 < a < 1, we have (1 + up < 1 + u for u > 0 so that 

|r(a; + l,x)| < / e““^(l + u) du < f-+ 

Jo \X X J 

< 2 e-^x"+^ (X > !)• 


When a > 1, we have 


poo 

|r(cu + l,x)| < e-^x“+'/ {e-pi + u)rdu 

Jo 

poo 

< / e““(l + u)du = 2e 

Jo 
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the first inequality holding when a < x and the second inequality when a > 1, since the 
integrand e““(l + u) < 1 on [0, oo). The resulting bound therefore holds for 1 < a < y. 
Collecting together these results, we therefore obtain the upper bound 

|r(w +l,x)| < 2e“^x“+^ (-l<a<x) (A.l) 


provided x > 1- 

Proceeding in a similar manner for the lower incomplete gamma function, we have 


|7(‘^ + 1,X)I = 






Jo 


< e 


-X^a+l 

Jo 


It is easily seen that 0 < e“^(l — n)“ < 1 on [0,1] when a > X) and hence we obtain the 
bound 

l7(w + l,x)l < (a>x)- (A. 2 ) 
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